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Abstract 

We propose a versatile Monte-Carlo method for pricing and hedg- 
ing options when the market is incomplete, for an arbitrary risk cri- 
terion (chosen here to be the expected shortfall), for a large class of 
stochastic processes, and in the presence of transaction costs. We il- 
lustrate the method on plain vanilla options when the price returns 
follow a Student-t distribution. We show that in the presence of fat- 
tails, our strategy allows to significantly reduce extreme risks, and 
generically leads to low Gamma hedging. Similarly, the inclusion of 
transaction costs reduces the Gamma of the optimal strategy. 

1 Introduction 

In their seminal 1973 article pQ, Black and Scholes (BS) have founded 
the very basis of modern financial mathematics. Their work has since 
been much studied and refined, and has become a rather abstract 
conceptual framework, deeply related with modern Probability The- 
ory [21 El • The BS model is the paradigm of complete markets, where 
every contingent claims can be replicated by a portfolio of underlying 
assets. Using a no-arbitrage principle one can deduce that any option 
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has a unique price, independent of the agent's risk preferences, which 
is given by the price of the replicating (or hedging) strategy. Mathe- 
matically speaking, these properties are equivalent to the existence of 
a unique equivalent martingale measure (also called the risk-neutral 
measure) under which one should average the final pay-off of an option 
to obtain its price. This measure is in general different from the 'true' 
(or objective) real world probability measure (under which we observe 
the evolution of financial assets) [21 IS]- The knowledge of this true 
probability distribution is thus in principle of no use for the pricing 
of options, although, as discussed in detail in jUEj, this message is in 
fact rather misleading. Within the BS framework, analytical formu- 
lae for the price and hedge exist for several cases, such as the simple 
European options. On the other hand, for more complex products 
like American or path dependent options, numerical procedures often 
have to be used [SJ [7] . 

However, the hypotheses of the BS model (Gaussian distribution 
and independence of log-returns, continuous time, absence of friction, 
etc.) have been widely questioned by practitioners and one can now 
observe a growing interest in the academic community for more general 
stochastic processes (Levy processes [5J [HI El E2] j stochastic volatility 
dUUaiinillllin], multifractal processes [El El HB]) which usually 
result in much more complex incomplete markets. In this case con- 
tingent claims cannot in general be replicated using underlying assets 
and the no-arbitrage principle is no longer enough to yield a unique 
price. One has to introduce additional criteria related to the risk pref- 
erences of the agents: the absence of perfect replication implies the 
existence of an intrinsic residual risk 0. 

In this framework one simple and popular methodology is the 
variance- minimizing hedging strategy |19l [3j. It consists in finding 
the (self-financing) portfolio whose difference with the pay-off of the 
option at maturity has the smallest variance. One drawback of this 
approach is that the risk-function is quadratic and therefore penalizes 
both profits and losses. Also, a quadratic measure of risk does not 
strongly penalize extreme risks. Alternative criteria based on higher 
moments of the distribution or on the notion of Value at Risk (VaR) 
and its extensions [201 have been recently considered in the literature 
|21l 1221 1231 [S] . In [5J the authors are interested in finding, given an 
initial investment, the portfolio strategy which maximizes the prob- 
ability of a successful hedge. Treated in a general setting and in a 
rather abstract way, this problem is shown to reduce to the replica- 
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tion of a particular knockout option. Although very appealing this 
solution can be hard to implement from a practical point of view. In- 
deed an explicit form of the option is only available in specific models 
and even in that case, as highlighted by the authors, the practical 
replication of such an option is not an easy task. The approach of 
\'2'2\ is much more concrete and numerically oriented. Nevertheless 
it is unsatisfactory from an optimization point of view: it is a static 
approach (the weights are determined at the initial time and remain 
constant) not well adapted to determine the full dynamical replication 
strategy. A similar observation can be made concerning the work pre- 
sented in |2S1, where the optimal static strategy that minimizes the 
fourth moment of the profit and loss distribution is determined. The 
approximate dynamical strategy is then constructed by 'translating' 
(in time) the optimal static strategy. An interesting observation made 
in that work is that the hedging strategy varies less rapidly with the 
underlying than the quadratic hedge, implying lower transaction costs 
(see also [3]). 

The aim of this paper is to propose a general Monte-Carlo algo- 
rithm that allows to price and hedge options using an arbitrary (but 
sufficiently smooth) risk criterion, and for a large class of stochastic 
processes describing the underlying. The algorithm, that generalizes 
the work of 4], is both easy to implement and versatile, and can be 
used for pricing different types of options. The paper is organized as 
follows: we first explain, following [I], the main ideas of our methodol- 
ogy. We then present some numerical results for different underlying 
processes, which show that extreme risks can be efficiently reduced 
compared to the standard BS hedge. As in |23j . we find that these 
extreme hedges have a smaller 'Gamma' at and around the money. 
Finally, we show how transaction costs can be treated within our 
method, and discuss several other possible extensions of our scheme. 

2 Description of the method 
2.1 Notations 

For simplicity we consider the case of a European option with one 
underlying asset of maturity T = Nr, where N is the number of re- 
balancing dates and r the time interval two dates. We denote the price 
of the underlying asset at time t& = kr hy x^, the strike by K and the 
final pay-off is [xn—K) + = max(xAr— K, 0). We suppose that the price 
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of the option only depends on the current price x k of the asset and 
call it by Ck(xk) at time tk- [If the volatility was stochastic, we should 
assume that the option price also depends on the current level of 
volatility a k and rather write Ck(xk,crk)]- The interest rate is assumed 
to be constant and equal to r and we define p = rr. Averaging 
(denoted by angled brackets (...)) will in the following always refer to 
the objective (real world) probability measure under which we observe 
the distribution of the asset returns, and not any abstract risk neutral 
measure. 

2.2 Principles 

The method we investigate here is an extension of the hedged Monte- 
Carlo strategy presented in 0J. Our aim is to construct a self-financing 
portfolio, whose wealth variation only depends on the variation of the 
asset price j3j, that best minimizes the chosen (instantaneous) risk 
measure. We denote by 4>k( x k) the fraction of the underlying asset in 
the portfolio at time k, when the asset price is x\~. Between time tk 
and tk+i the self-financing condition leads to a local wealth balance 
given by: [3] 

AW k = e p C k (x k ) - C k+l (x k+ i) + (j) k (xk)(xk+i - e p x k ) (1) 

The measure of the quality of the replication is given by a local risk 
function U{AW k )- The average risk, over all paths of the real process, 
is thus given by: 



For purposes of illustration, we have chosen in the following a function 
U(AWk) which penalizes losses that exceed a certain threshold — Aq: 



where the exponent q was chosen to q = 1, corresponding to what 
is called an unconditional expected shortfall. The generalization to 
arbitrary q, or in fact to other functional forms for U(AWk), does not 
lead to any numerical difficulty (q > 1 penalizes more strongly extreme 
losses). The minimization of IZk is quite sensible from a financial point 
of view: one tries to control in a marked to market way, during the 
whole life of the option, the occurrence of downside moves. Choosing 
a large negative Aq means that one aims at controlling extreme losses. 




(2) 



U(AW k ) = (A - AW k )\ 



A - AW k \ q l AWk</ll 



■0 ' 



(3) 
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From a practical point of view we solve the above optimization 
problem using Monte Carlo simulations, that allows one to use a rather 
general stochastic process for the price evolution. We generate Nmc 
trajectories of the asset price over which we will average. Following 
PUIH> we decompose the functions Ck(x) and 4>k(x) on a set of p fixed 
basis functions: 



Ck(x) = £tJC*(*) (4) 
a=l 

M*) = (5) 



a=l 

Doing this we reduce the original functional optimization (find the 
functions (pk and Cj.) to a numerical optimization: we now have a 
minimization problem in term of the parameters 7^ and ip\. If p is 
large enough we expect to have a good approximation of the true 
functional solution (see the following section for numerical implemen- 
tation). We solve the problem by working backward in time from 
maturity, where the option is worth its known final pay-off. For each 
time k, we decompose the problem into the following steps: 

• If the time discretisation mesh is sufficiently small, one can ap- 
proximate Ck(xk) by Ck+i{xk) whose functional form is already 
known from the previous step. 1 We then find the coefficients ip\ 
which minimize the average risk over the Nmc paths: 

Nmc , p \ 

n = E A o-(e^ +1 (4)-c fe+1 (4 + o)-(4+i-^4)E^ i?, a(4) 

t=\ ^ a=l ' 

This can be done using a steepest gradient method. Indeed we 
can easily compute the partial derivative of lZ* k with respect to 
the coefficients tp k a : 

&R* Nmc 

■q-t; = ~ E ( X k+1 ~ ePx k) F a ( a:; fc) 1 AW^<Ao- 

• Using the fact that, on average, the local wealth balance must 
equal zero, we now compute the coefficients 7^ by solving the 



1 A better approximation that takes into account the (already known) time derivative 
of C k (x k ) is to write C k {x k ) ~ 2C k +i{x k ) ~C k+2 (x k ). 
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least square problem: 



min 

7 



E 

i=i 



E 

.0=1 



■r i 



Ch 



+H x fc+1/ 



M x i)( x k+i 



which is easily done using standard procedures |25j . 



3 Numerical results 
3.1 Implementation issues 

In this section we compare the results obtained following a standard 
Black and Scholes strategy (delta hedging) with those obtained follow- 
ing strategies with different values of the threshold Aq, as explained 
above. We price a European option, with a maturity of 1 year and 
an annualized volatility a = 20%. We choose a rather small number 
of time intervals when re-hedging is possible, N = 10. The initial 
stock price is xo = 100. We use Nmc = 20000 trajectories for averag- 
ing. We first consider realizations of a standard geometric Brownian 
motion with a constant drift /x = 0.05: 

dx t = x t (fJ-dt + adWt) 

It is well known that financial time series are very poorly represented 
by such a process and display much heavier tails. To qualitatively 
account for this fact we also use realizations of a fat tailed process 
where we replaced the previous Brownian motion Wt by a fat tailed 
process L t whose increments are distributed according a Student— t 
distribution, with v = 4 or v = 6 degrees of freedom: 2 

dx t = Xt (fJ.dt + crdL t ) 

The value of v characterizes the power-law decay of the distribution 
for large arguments; v = 4 is in the range of reported value for this 
exponent for rather liquid markets. The value v = 6 corresponds 
to faster decaying tails; the limit v — > oo corresponds to the usual 
Black-Scholes model. Some markets (like emerging country markets, 

See | 26 | for option pricing with Student— t distribution of returns on all time scales, 
to be contrasted with the present model where only the shortest time scale returns are 
distributed according to a Student— t distribution. 
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or emerging country currencies) would correspond to small values of v 
(for example v ~ 1.5 for the Mexican Peso). We use p = 20 basis func- 
tions, which we find to be accurate enough. Following 4, we choose for 
F^ piecewise linear functions and for C% piecewise quadratic functions 
with both the same adaptive breakpoints. These breakpoints are cho- 
sen so that at each stage the same number of simulated trajectories, 
Nmc/{p + 2), fall between two successive breakpoints. is worth 
below the a th breakpoint, 1 above the a + 1 th breakpoint and is linear 
between these two values. C\ is taken as the integral of Fjf which is 
worth below the a th breakpoint. Our numerical simulations were 
systematically conducted as follows: we first generate a set of trajec- 
tories and apply our algorithm to find the coefficients 7^ and ip%, i.e., 
the price of the option and the optimal hedge. We then simulate a 
new set of paths to compute different statistical indicators of the per- 
formance of the proposed strategies. In other words, different paths 
are used for the optimization, and for backtesting the optimization in 
an 'out of sample' fashion. 

3.2 Expected shortfall hedging in the Black- 
Scholes case 

We first present the results obtained within the framework of the Black 
and Scholes model. Since perfect replication is theoretically possible 
in this case, we expect that the minimum short-fall strategy we find 
should be close to the Black and Scholes strategy. In the continuous 
time limit, the Black-Scholes strategy actually leads to a zero expected 
shortfall, for any value of Ao < 0. This is indeed what we observe in 
Fig. ^where we plot, for a strike price equal to K = 110, the optimal 
solutions found with different values of Ao and the BS strategy: they 
all look very similar, in particular for small values of Ao (—0.5 and 

-I)- 

This observation is confirmed by the shape of the distribution of 
the final wealth at the maturity date of the option, for the different 
strategies (Fig. |2J). 

3.3 The case of a fat-tailed dynamics 

We are now interested in a market whose dynamics is governed by a fat 
tailed Levy process, where the relative price increments are indepen- 
dent identically distributed with jumps. In this case the market is no 
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Figure 1: Optimal number of risky assets in the hedging portfolio, as a 
function of the level x of the underlying asset for different strategies within 
the Black-Scholes. As expected the different curves are very similar. The 
value of the strike is K = 110, and the time is half the maturity of the 
option: k = N/2. 
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Figure 2: Distribution of the final wealth for different strategies for a BS 
market (upper left: BS, upper right: A = 0, lower left: A = —0.5, lower 
right: Ao = —1)- The four curves are not very different from one another. 
The value of the strike is K = 110. 
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longer complete (existence of unhedgeable jumps) and perfect replica- 
tion does not exist. The use of a subjective criteria is then needed for 
pricing and hedging purposes, we now expect to get different strate- 
gies depending on the value of Aq. Our expectations are numerically 
confirmed in Fig. |3] where we clearly see strong differences between 
the strategies. In particular, we observe that extreme losses hedging 
(corresponding to a large value of |Ao|) leads to a flatter function 4>(x) 
(i.e. a smaller Gamma). This was already emphasized in \2'A[ 15] and 
can be quite interesting in presence of transaction costs (see section 
4). 




Figure 3: Optimal number of risky assets in the hedging portfolio, as a 
function of the level x of the underlying asset for different strategies in the 
case of a fat tailed market (left: v = 6 and right: v = 4). The hedging 
strategies now present a very different dependence on the underlying, in 
particular when |Ao| increases. The value of the strike is K — 110 and the 
time is half the maturity of the option: k = N/2. 

Using an independent set of paths, we can check that our method 
indeed leads to smaller local expected shortfalls when the correspond- 
ing optimal strategy is adopted, as can be seen in Tables Q an d EJ 
where we give both the unconditional expected shortfall 1Z, and the 
conditional expected shortfall (noted ESF), defined asTZ/V, where V 
is the probability to exceed the threshold Ao- These expected short- 
falls are computed between two re-hedging times k and k + 1, where 
k = N/2 corresponds to half the life of the option. 
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strategy 


BS 


A = 


A = -5 


A = -10 


ESF(O) 


-0.77 


-0.88 


-0.93 


-1.19 


ESF(-5) 


-2.26 


-2.32 


-2.17 


-1.98 


ESF(-IO) 


-4.38 


-4.00 


-3.05 


-2.21 


K(0) 


-0.22 


-0.25 


-0.33 


-0.44 


TZ(-5) 


-0.011 


-0.014 


-0.010 


-0.015 


K(-IO) 


-0.002 


-0.003 


-0.002 


-0.001 



Table 1: Conditional expected shortfall (ESF) and unconditional expected short- 
fall 1Z, for different threshold Ao and different strategies, in the case v = 6. These 
local quantities are computed for a time equal to half the maturity of the option 
(here 1 year), for a strike K = 110. 



strategy 


BS 


A = 


A = -5 


A = -10 


ESF(0) 


-0.88 


-1.04 


-1.04 


-1.10 


ESF(-5) 


-4.45 


-4.40 


-3.94 


-3.15 


ESF(-IO) 


-8.58 


-8.60 


-7.76 


-5.90 


71(0) 


-0.23 


-0.25 


-0.35 


-0.38 


n(-5) 


-0.031 


-0.035 


-0.028 


-0.028 


n(-io) 


-0.015 


-0.016 


-0.011 


-0.009 



Table 2: Conditional expected shortfall (ESF) and unconditional expected short- 
fall 1Z, for different threshold Ao and different strategies, in the case v = 4. These 
local quantities are computed for a time equal to half the maturity of the option 
(here 1 year), for a strike K = 110. 

We can also check that our method leads to satisfactory results for 
global quantities (i.e. concerning the wealth balance at the end of the 
option lifetime). We therefore determined the distribution of the final 
wealth (Fig. 0] and . We clearly see that the strategy proposed here 
can significantly reduce the value of the extreme losses (note that 
because of the power-tails of the return distribution, these extreme 
losses can still be large). Moreover we observe a remarkable change in 
the shape of this distribution: from relatively peaked for small |Ao| 
but with an appreciable number of extreme events, it becomes broader 
but with a smaller support when |Ao| increases. This fact can be 
qualitatively explained: when |Aq| is large, the constraint is easier to 
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fulfill but losses of amplitude less than |An| are not penalized, leading 
to a broader looking but more sharply truncated final distribution. 




-50 -40 -30 -20 -10 10 -50 -40 -30 -20 -10 10 

600 | 1 1 1 1 1 1 450 I 1 1 1 1 1 1 




Figure 4: Distribution of the final wealth for different strategies in the case 
of a fat-tailed market with v — 6 (upper left: BS, upper right: Ao = 0, lower 
left: A = -5, lower right: A = -10). 

We give in Tables 01 and |I] the mean and the standard deviation 
of these final wealth distributions, as well as the initial price of the 
option and the associated value at risk and expected shortfalls. Notice 
that the option price is smaller than the Black-Scholes price, which is 
expected when the moneyness is small: non zero kurtosis indeed leads 
to a decrease of the at-the-money volatility [3] as can be seen in Fig. 
El The option price in fact decreases when larger risks are hedged. 
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Figure 5: distribution of the final wealth for different strategies in the case 
of a fat-tailed market with v — 4 (upper left: BS, upper right: A = 0, lower 
left: A = —5, lower right: A = —10). 
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Figure 6: implied volatility smile at t — 0. v — 4 and So = 100 
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strategy 


BS 


A = 


An = -5 


A = -10 


rvntion nrirp 


5.29 


5.15 


5.14 


5.08 


mean of final wealth 


0.15 


0.01 


0.03 


0.04 


std of final wealth 


3.11 


3.23 


3.41 


4.01 


VaR 0.1% 


-21.81 


-20.77 


-19.04 


-18.40 


ESF 0.1% 


-7.96 


-8.50 


-7.63 


-6.73 


VaR 1% 


-10.04 


-10.63 


-8.82 


-9.47 


ESF 1% 


-4.88 


-4.91 


-4.41 


-3.69 


VaR 5% 


-4.85 


-5.39 


-5.28 


-6.19 


ESF 5% 


-3.39 


-3.46 


-2.56 


-2.28 



Table 3: Statistical characteristics of the global wealth balance for v — 6. 
The strike is K = 110. 



strategy 


BS 


A = 


A = -5 


A - -10 


option price 


5.29 


5.08 


4.99 


4.89 


mean of final wealth 


0.29 


0.09 


0.07 


0.09 


std of final wealth 


4.02 


4.12 


4.20 


4.82 


VaR 0.1% 


-34.90 


-36.49 


-28.37 


-27.45 


ESF 0.1% 


-16.14 


-15.21 


-14.05 


-9.87 


VaR 1% 


-13.59 


-14.28 


-12.08 


-12.72 


ESF 1% 


-9.34 


-9.34 


-7.85 


-6.06 


VaR 5% 


-5.58 


-6.19 


-6.32 


-7.76 


ESF 5% 


-5.48 


-5.55 


-4.23 


-3.51 



Table 4: Statistical characteristics of the global wealth balance for v — 4. 
The strike is K = 110. 

4 Application to transaction costs 

As mentioned above, hedging against extreme risks generically leads 
to a strategy that varies more slowly with the underlying asset price. 
This can be of great interest in the presence of transaction costs. 
These costs can in fact be endogenously taken into account within 
the present numerical scheme, which allows to determine how both 
the price and the optimal hedge are impacted by transaction costs. 
Previous analytical work on this problem in the framework of the BS 
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model can be found in |27l 17). and further discussions and extensions 
to the case of non Gaussian markets can be found in |28| . 

We model friction by adding to the wealth balance a cost propor- 
tional to the number of bought or sold assets and to its price. EqQ 
now becomes 

A W k = e p C k (x k ) - Cfc+i (xfc+i) + 4> k (x k ) (x k +i - e p x k ) - (3x k \ 4>k - <t>k-x\ , 

(6) 

where (3xk represents the transaction costs per share. Following the 
same steps as above, we now want to minimize the risk function 

Nmc / 

IV k = E A ° " { ePC k+i(4) - Cfc+i(4+i)) - - e p xl)M* 

o 1 1 I I i / I \ 
+P x k\ x k ~ X k-l\-Q^\ X k) 

where we approximated \4>k — <fik-i\ by ^|Ax|. This is justified if 
the time step is sufficiently small, and is the key step to make the 
problem tractable. However, since this involves the derivative of 
we have preferred to work with a smooth parameterization of the 
function <pk with only two optimization parameters, rather than the 
full decomposition over a set of basis functions, as was used above. 
We have checked that the following choice gives very similar results 
than the ones obtained above in the absence of transaction costs. We 
thus take: 

(t) k {x) = - (l + [tanh|^ fc 7W fc |] A .sign(7W fc )) 

with a rescaled moneyness Aik given by: 

x - Ke-^-^ 

M k = 



a^T - t k 

By varying the two variational parameters Ak and (3k, we can, at 
each time step, optimize any risk measure. Using essentially the same 
numerical optimization procedure as above we are able to find the 
parameters Ak and (3k, and thus the optimal strategy (j)* k . In order 
to obtain the option price, we have then to solve the following least 
square problem: 
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mm 

7 



a=l ^ 



e p 4) 



+P x k\ x k ~ x t-l\~Qz(. x k. 



We have numerically tested the above scheme in the case of a BS 
market with the same characteristics as in the previous section, and 
for different values of the friction parameter j3. We compared our 
results with those obtained following a naive BS strategy and the more 
advanced Leland strategy. Using simple arguments, similar in spirit 
to the above approximation on A<p, Leland showed in [27 how the 
BS strategy can be modified to account for transaction costs. Indeed 
using a modified volatility 



2 1 

a L = a\ 1 + 2/3 



7T aV At 

instead of the real volatility a leads to a strategy which, on average, 
approximately covers the transaction costs and hedges the risk. Since 
<7£ > cr, the option price is, as expected, higher than in the BS case. 
Figure [7| compares the obtained optimal strategies using our method 
with both the BS and Leland hedging schemes. As expected, the costs 
affects the BS strategy in such a way as to reduce its at the money 
Gamma. Tables El El and El give a summary statistics of our results. 
We note that even in the case where the stock price is log-normal, our 
strategy allows to improve significantly over the Leland strategy if the 
threshold |Ao| is large enough: compare Tables El and El Using the 
optimal strategy allows one to simultaneously reduce the occurrence 
of large risks (measures both by the VaR and the ESF) while keeping 
the option price lower than in the Leland scheme. 



5 Conclusion 

In this paper we have extended the work of [I] and proposed a general 
numerical (Monte-Carlo) methodology for the pricing and hedging of 
options when the market is incomplete, for an arbitrary risk crite- 
rion (chosen here to be the expected shortfall) and in the presence 
of transaction costs. We have shown that in the presence of fat-tails, 
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50 75 100 125 150 175 

x 

Figure 7: Optimal number of risky assets in the hedging portfolio, as a 
function of the level x of the underlying asset for different strategies: Black- 
Scholes, Leland, A = —1 and A = —5, for (3 = 0.05. The value of the 
strike is K=110. 
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p 


0.005 


0.01 


0.05 


nntinn nripp 

\JLJ\JL\JL± I V 1 IV. V. 


5.29 


5.29 


5.29 


mean of final wealth 


-0.41 


-0.90 


-4.78 


std of final wealth 


2.38 


2.45 


3.28 


VaR 0.1% 


-10.95 


-11.86 


-19.91 


ESF 0.1% 


-1.45 


-1.47 


-1.33 


VaR 1% 


-7.32 


-8.06 


-14.83 


ESF 1% 


-1.57 


-1.67 


-2.07 


VaR 5% 


-4.47 


-5.15 


-11.14 


ESF 5% 


-1.75 


-1.84 


-2.25 



Table 5: Impact of the transaction costs for the B&S strategy. 



p 


0.005 


0.01 


0.05 


option price 


5.69 


6.08 


9.27 


mean of final wealth 


0.10 


0.13 


0.32 


std of final wealth 


2.39 


2.44 


2.98 


VaR 0.1% 


-9.83 


-9.79 


-10.32 


ESF 0.1% 


-1.54 


-1.53 


-1.28 


VaR 1% 


-6.47 


-6.47 


-6.83 


ESF 1% 


-1.55 


-1.54 


-1.60 


VaR 5% 


-3.80 


-3.84 


-4.27 


ESF 5% 


-1.67 


-1.64 


-1.62 



Table 6: Impact of the transaction costs for the Leland strategy. 

our strategy allows to significantly reduce extreme risks, and gener- 
ically leads to low Gamma hedging, as anticipated in |23l 15]. Many 
other risk criteria could be considered, in particular functions that 
give more weights to extreme losses. We focused in this work on plain 
vanilla European options, but (as shown in [I]) the method is readily 
extended to a large family of exotic options. Finally, we showed how 
our method allows to deal consistently with transaction costs. When 
compared to the standard Leland hedging scheme, our optimal strat- 
egy leads both to lower option prices and better hedging of large risks, 
even in the simplest case of a log- normal market. 
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p 


0.005 


0.01 


0.05 


nntinn nripp 

\JLJ\JL\JL± I V 1 IV. V. 


5.77 


6.22 


8.30 


mean of final wealth 


0.00 


-0.04 


-0.26 


std of final wealth 


2.39 


2.44 


3.09 


VaR 0.1% 


-11.05 


-10.95 


-10.87 


ESF 0.1% 


-1.34 


-1.45 


-1.27 


VaR 1% 


-7.05 


-7.13 


-7.54 


ESF 1% 


-1.61 


-1.54 


-1.45 


VaR 5% 


-4.01 


-4.11 


-5.09 


ESF 5% 


-1.84 


-1.80 


-1.51 



Table 7: Impact of the transaction costs using the optimal strategy, where 
these costs are accounted for. A = — 1 



p 


0.005 


0.01 


0.05 


option price 


5.63 


5.99 


8.44 


mean of final wealth 


0.00 


-0.05 


-0.75 


std of final wealth 


2.66 


2.74 


3.39 


VaR 0.1% 


-8.92 


-9.17 


-11.95 


ESF 0.1% 


-1.17 


-1.27 


-1.44 


VaR 1% 


-6.12 


-6.37 


-8.60 


ESF 1% 


-1.28 


-1.31 


-1.48 


VaR 5% 


-4.06 


-4.25 


-6.07 


ESF 5% 


-1.30 


-1.33 


-1.56 



Table 8: Impact of the transaction costs using the optimal strategy, where 
these costs are accounted for. A = —5 

There are many extensions of the above method that would be 
worth investigating, in particular the case where the underlying has a 
stochastic volatility with some persistence, such as, for example, the 
models studied in [TJ1 [EH El EEl HE]- In this case, both the price 
and the optimal hedge should explicitly depend on the local value of 
the volatility, or of a noisy estimate of this volatility. Other hedging 
instruments, like options of different maturities, could in this case be 
included in the local wealth balance to reduce the risk further. An- 
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other interesting path is to consider the problem of hedging a whole 
portfolio of options. Contrarily to the Black-Scholes case, where the 
optimal strategy for the whole portfolio is the linear sum of the indi- 
vidual hedges, extreme value hedges lead to a non linear composition 
of the individual hedges. 
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